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Parton distributions in the virtual photon target are investigated in perturbative QCD up to 
the next-to-next-to-leading order (NNLO). In the case A 2 <C P 2 <C Q 2 , where — Q 2 (-P 2 ) is 
the mass squared of the probe (target) photon, parton distributions can be predicted completely 
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up to the NNLO, but they are factorisation-scheme-dependent. We analyse parton distributions 
in two different factorisation schemes, namely MS and DIS 7 schemes, and discuss their scheme 
dependence. We show that the factorisation-scheme dependence is characterised by the large-x 
behaviours of quark distributions. Gluon distribution is predicted to be very small in absolute 
value except in the small-x region. 
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FIG. 1: Deep inelastic scattering on a virtual photon in the e + e collider experiments. 

I. INTRODUCTION 

It is well known that, in e + e" collision experiments, the cross section for the two-photon 
processes e + e~ — > e + e~ + hadrons illustrated in Fig.l dominates at high energies over other 
processes such as the one-photon annihilation process e + e~ — > 7* — > hadrons. The two- 
photon processes are observed in the double-tag events where both the outgoing e + and e~ 
are detected. Especially, the case in which one of the virtual photon is far off-shell (large 
Q 2 = — q 2 ), while the other is close to the mass-shell (small P 2 = —p 2 ), can be viewed as 
a deep-inelastic scattering where the target is a photon rather than a nucleon [1]. In this 
deep-inelastic scattering off photon targets, we can study the photon structure functions, 
which are the analogues of the nucleon structure functions. 

A unique and interesting feature of the photon structure functions is that, in contrast 
with the nucleon case, the target mass squared P 2 is not fixed but can take various values 
and that the structure functions show different behaviours depending on the values of P 2 . 

In the case of a real photon target (P 2 = 0), unpolarised (spin-averaged) photon structure 
functions F^XjQ 2 ) and F2(x,Q 2 ) were studied first in the parton model [2], and then 
investigated in perturbative QCD (pQCD). In the framework based on the operator product 
expansion (OPE) [3] supplemented by the renormalisation (RG) group method, Witten [4] 
obtained the leading order (LO) QCD contributions to F 2 7 and and, shortly after, the 
next-to-leading order (NLO) QCD contributions were calculated by Bardeen and Buras [5]. 
The same results were rederived by the QCD improved parton model powered by the parton 
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evolution equations [6, 7]. Recently, some of the parameters which are necessary to evaluate 
the next-to-next-to-leading order (NNLO) corrections to F 2 7 were calculated and also the 
NNLO solution of the evolution equations for the real photon case was dealt with [8]. When 
polarised beams are used in e + e~ collision experiments, we can get information on the spin 
structure of the photon. The QCD analysis of the polarised photon structure function 
gJ(x,Q 2 ) for the real photon target was performed in the LO [9] and in the NLO [10, 11]. 
For more information on the recent theoretical and experimental investigation of unpolarised 
and polarised photon structure, see the review articles [12]. 

For a virtual photon target (P 2 ^ 0), we obtain the virtual photon structure functions 
F^(x,Q 2 ,P 2 ) and F2(x,Q 2 , P 2 ). They were studied in the LO [13] and in the NLO [14] by 
pQCD. In fact, these structure functions were analysed in the kinematical region, 

A 2 <P 2 <Q 2 , (1.1) 

where A is the QCD scale parameter. The advantage of studying a virtual photon target 
in this kinematical region (1.1) is that we can calculate the whole structure function, its 
shape and magnitude, by the perturbative method. This is contrasted with the case of the 
real photon target where in the NLO and beyond there appear nonperturbative pieces. The 
virtual photon structure functions and Fl were also studied by the QCD improved parton 
model [15-18]. In the same kinematical region (1.1), the polarised virtual photon structure 
function gJ(x,Q 2 , P 2 ) was investigated up to the NLO in pQCD [19], and its target mass 
effects were studied [20]. Moreover, the first moment of g1(x,Q 2 , P 2 ) was calculated up to 
the NNLO [21]. 

Recently the three-loop calculations were made for the anomalous dimensions of both 
quark and gluon operators [22, 23] and also for the photon-quark and photon-gluon split- 
ting functions [24]. Using these results, we analysed the virtual photon structure functions 
F^(x, Q 2 , P 2 ) and F2(x, Q 2 , P 2 ) in the kinematical region (1.1) in QCD. We could give defi- 
nite predictions for P 2 7 (x, Q 2 , P 2 ) up to the NNLO and for F2(x, Q 2 , P 2 ) up to the NLO [25], 
and also for the target mass effects on these functions [26]. 

In parton picture, structure functions are expressed as convolutions of coefficient functions 
and parton distributions in the target. And these parton distributions will be used for 
predicting the behaviours of other processes. When the target is a virtual photon with P 2 
being in the kinematical region (1.1), then a definite prediction can be made for its parton 
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distributions in pQCD. However, it is well known that parton distributions are dependent on 
the scheme which is employed to factorise structure functions into coefficient functions and 
parton distributions. It is possible that parton distributions obtained in one factorisation 
scheme may be more appropriate to use than those acquired in other schemes. Indeed, 
for the case of the polarised virtual photon target, the NLO QCD analysis was made for 
its parton distributions in several factorisation schemes and their scheme dependence was 
closely studied [27]. 

In this paper we perform the QCD analysis of the parton distributions in the unpolarised 
virtual photon target up to the NNLO 1 . We investigate the distributions of the flavour 
singlet and nonsinglet quarks and gluons in two factorisation schemes, namely MS and DIS 7 
schemes 2 . Using the framework of the QCD improved parton model, we give, in the next 
section, the explicit expressions for the moments of the flavour singlet (nonsinglet) quark and 
gluon distributions up to the NNLO. In Sec. Ill, we show first that the parton distributions 
obtained up to the NNLO are independent of the renormalisation-prescription which is 
chosen in defining the QCD coupling constant a s . And then we explain the two factorisation 
schemes, MS and DIS 7 schemes, which we consider in this paper. The behaviours of the 
parton distributions near x — 1 and their factorisation-scheme dependence are discussed 
in Sec. IV. The numerical analysis of the parton distributions predicted by MS and DIS 7 
schemes will be given in Sec.V. The final section is devoted to the conclusion. 

II. PARTON DISTRIBUTIONS IN THE VIRTUAL PHOTON 

In this section the parton distributions in the virtual photon are analysed in the framework 
of the QCD improved parton model [29] with the DGLAP parton evolution equations. 

Let q±(x, Q 2 , P 2 ), G±(x, Q 2 , P 2 ), T±(x, Q 2 , P 2 ) be quark with i-flavour, gluon, and pho- 
ton distribution functions with ± helicities in the virtual photon with mass — P 2 . Then the 

1 This work is a sister version of Ref . [27] . 

2 Part of the results in this paper was reported at the 8th International Symposium on Radiative Corrections 
(RADCOR)[28]. 
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spin-averaged parton distributions in the virtual photon are defined as 



(2.1) 



In the leading order of the electromagnetic coupling constant, a = e 2 /An, T 7 does not evolve 
with Q 2 and is set to be 

V(x,Cf,P 2 ) =S(l-x) . (2.2) 

Instead of using the quark distributions q\ it is advantageous to treat the flavour singlet 
and nonsinglet combinations defined as follows: 



?2 = £ 9* , 



(2.3) 
(2.4) 



where n/ is the number of relevant active quark (i.e., the massless quark) flavours and e$ is 
the electromagnetic charge of the active quark with flavour % in the unit of proton charge. 
Now introducing a row vector 



qi = G\ q NS ), 



(2.5) 



the parton distributions q^(x, Q 2 , P 2 ) in the virtual photon satisfy an inhomogeneous evo- 
lution equation [6-8] 



dq^ x ,Q 2 ,p 2 ) = kM 2 )+ ! ld y<? M 2 ,P 2 )p£,Q 2 ) , 

Jx y y 



(2.6) 



d InQ 2 

where the elements of a row vector k = (ks, k G , k^s) refer to the splitting functions of 7 to 
the singlet quark combination, to gluon and to nonsinglet quark combination, respectively. 
The 3x3 matrix P(z, Q 2 ) is expressed as 



P(z,Q 2 ) 



P q s q (z,Q 2 ) P Gq (z,Q 2 ) 



\ 



(2-7) 





P qG (z,Q 2 ) P GG (z,Q 2 ) 
V P q N q s (z,Q 2 )J 

where Pab is a splitting function of 5-parton to A-parton. 

There are two methods to obtain the parton distribution q y (x, Q 2 , P 2 ) up to the NNLO. 
In one method, we use the parton splitting functions up to the three-loop level and we solve 



numerically q y (x, Q 2 , P 2 ) in Eq.(2.6) by iteration, starting from the initial quark and gluon 
distributions of the virtual photon at Q 2 = P 2 . The interesting point of studying the virtual 
photon with mass — P 2 is that when P 2 3> A 2 , the initial parton distributions of the photon 
are completely known up to the two-loop level in QCD. The other method, which is more 
common than the former, is by making use of the inverse Mellin transformation. From now 
on we follow the latter method. First we take the Mellin moments of Eq.(2.6), 

dq-y(n,Q 2 ,P 2 ) 



d k(n,Q 2 ) • q(n.(f.r 2 )r(n.(f) , 

where we have defined the moments of an arbitrary function f(x) as 



(2.8) 



f(n)= r 

Jo 



dxx^fix). 



(2.9) 



Henceforce we omit the obvious n-dependence for simplicity. Expansions are made for 
the splitting functions k(Q 2 ) and P(Q 2 ) in powers of the QCD and QED coupling constants 
as follows: 



a s (Q 2 



2vr 



P(Q 2 ) = ^flp^ + 



2tt 



2\1 2 



p(D + 



2tt 



(2.10) 
(2.11) 



and a new variable t is introduced as the evolution variable instead of Q 2 [30], 



f 4ln^ 



A> « S (Q 2 ) 



(2.12) 



The solution q' Y (t)(= q 1 (n, Q 2 , P 2 )) of Eq.(2.8) is decomposed in the following form: 



qi(t) = qiW(t) + q« 1 )(t) + qM(t) , 



7(1), 



iT( 2 )f 



(2.13) 



where the first, second and third terms represent the solution in the LO, NLO and NNLO, 
respectively. Then they satisfy the following three differential equations: 

dqi<®(t) 



a 



-kW + q«°\t)pW , 



dt a s (t) 



dt 



2ir 



dqi{ y = Rk{2) + r^i) 2 ^(o )(0 p 

dt 2tt 2tt K(2) V 2tt / y W ni) 

+ q 7(D (t) jRp(1) + q7(2) (t) p(0) ) 



(2.14) 
(2.15) 

(2.16) 



where 



■Kff(i) 


- fcw 

— rv 


Rp(i) 


= pt 1 ) 


Rk(2) 


= 


Rp(2) 


= P( 2 ) 



A, 

2(3o 

A 

2/V 



(2.17) 
(2.18) 



2/3 4 VVo^ /V 

where we have used the fact the QCD effective coupling constant a s (Q 2 ) satisfies 



(2.20) 



<7(P) I I 7(P)>U=P- = + ^4P^ (2) }, < = ^ (or ^),G,iV5 . (2.22) 



~dh^~ ~ - % ^r ~ Pl i^w~ " A i4^" + " ' ' (2 - 21) 

where (3q = 11 — |n/ and /3i and $2 were known [31, 32]. Note that the P 2 dependence of 
cf solely comes from the initial condition (or boundary condition) as we will see below. 

The initial conditions for q 1 ^, q 1 ^ and q 1 ^ are obtained as follows: For — p 2 = P 2 ^> 
A 2 , the photon matrix elements of the hadronic operators O l n (i = S (or ip), G, NS) can be 
calculated perturbatively. Renormalising at /i 2 = P 2 , we obtain at two- loop level 

The An 1 ^ and A 1 ^ 1 terms represent the operator mixing between the hadronic operators and 
photon operators in the NLO and NNLO, respectively, and the operator mixing implies that 
there exist parton distributions in the photon. Thus we have, at /i 2 = P 2 (or at t — 0), 

9*°>(0) = 0, ^)(0) = £A«, ^ (2) (0) = ^^Ai 2 ), (2.23) 

with 

A^=(A s J l \A^\A^) , 1 = 1,2. (2.24) 

Actually, the initial quark distributions emerge in the NLO (the order a) and gluon distri- 
bution in the NNLO (the order aa s ). The expressions of and A$ in MS scheme are 
enumerated in Appendix A 1. 
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With these initial conditions, the solutions q^ 0) {t), <? 7(1) (t) and q^ 2 \t) are given by 



An 

a s (t) 





a s (t) 







2P(") 
00 



(2.25) 



3 (0) 



a 



P(o) i 1 


~a s (t)~ 


0o | 




_a s (0)_ 





— Rk(i) + 2a Rp(i) 



-2 a \ I dre( p(0) -^ R e -P^r \ + q7 (i) (0) 



a,(0) 



00 



^2.26) 



« 7(2) (t) = 



2tt 



2^ Rk(2) + 2<2 -Rp(2) 



& + p(°) 



Q s (t) 
« s (0) 



=(0) 



A) 



2a| / dre^"^ i2 P(2) 



+ 



a 
27 



-Ri^(i) + 2a -Rp(i) 



p(o) 





1 - e p,,v 







i2 P(1) e-^ +p( °^ \e^ +pi0) » 

p(0), 



-2a|^ rfr £dr'e^ 0) -^y R P{1) e 

xe p(0) ^p (1) e -^ +p(0) ^}e^^ (0) )* 

2P(°) 

a s (t) 



+9 7(2, (0) 



a s (0) 



where 



2tt/3 1 - 

P0 



(2.27) 



(2.28) 



For the case of the real photon target, the NLO solution q 7<1 )(t) was given in [7] (see also 
[33]) and the NNLO solution q^ 2 \t) in [8]. 

The moments of the splitting functions are related to the anomalous dimensions of oper- 
ators as follows [25]: 



p(o) = _i^(o) p(i) = P (2) = 

A ^ In i A 8 16 

fc (0) = \k*\ fc« = fe< 2 > = 1k« . 



(2.29) 



The evaluation of q^°\t), q^(t) and q^ 2 \t) in Eqs.(2.25)-(2.27) can be easily done by 
introducing the projection operators P™ [5] 



(2.30) 



,NS 



pn pn 

i j 



E^ = 1 > 



(2.31) 



for % ^ j, 
P™ for % = j, 

where A™ are the eigenvalues of the matrix \ Then, rewriting a s (0) and a s (t) as a s {P 2 ) 
and a s (Q 2 ), respectively, we obtain 



a 



8tt(3 } a s (Q 2 ) 



l7r K (0)Ypn^_{ 1 _ 



«,(P 2 ) 



(2.32) 



q w/ L87taJ \ ™ 4" ^ A) n v ^ ^ 



pn ^C 1 ) pn i 
7 In r { 1 



pn ^(1) pr 
i /n 7 



X { 1 - 



a s (Q 2 ) 

MP 2 ) 



+ <^ ft-' ' V - J ' _ ^l|f(o) V P n 1 



x 1- 



a s (Q 2 ) 
a s (P 2 ) 



(2.33) 
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(*)/ 



a 



8np l 



■otsiQr 

An 



n W ^ 1 V 2)^ n p ^ 1 1 - df\ 2 J 



n A, 



pn ^1) pn X _ d n pn jg) pn j _ g + g 

^ 2/3 + A? - A™ 1 - + 4^ 4/3 + A™ - A™ 1 - 



pre ij( 2 ) pn 

+^ 0) E J ' 



-? 4/5 + A"-A™ l-< 



P" 7i 1} PJ! ^n> P % 



(1) 



1 



(2/?o-Ar + A^)(4/5 + A«-A«) 1-d? 



.ifW^v pn , K {i) v 5! )p ? 1 _ r-(2) v p n 1 

n A>4^ 1 n 4^ 2/3 + A"-A™ 1-d? n 4^ 1 



-(1) nn 



+ 2 ^" E 2 ^-^ - 2 ^'§ E t<c - 2 ^i 2) E t 



3,i 



00 

A, 



a s (Q' 



a s (P< 



d?-l 



+ J JffOlf^V pn ( 1 _ d n ) _ K (0)P] L y- P i^ n P 



A? < 



fty jyjiPP[_ (0) v g 7n X) iy 7^ 1 

A, 4^ 2/3 + A™ - A™ * £-< (2ft + A™ - A£)(2/? + A™ - A") df 



3,1 



n A> V " 4^ 2 Ao + A? - A™ < 



pn J^ 1 ) pn 
- r j 7n -fj 



-2AoA« e 



pn ^7(1) pn 
r t 7n -Tj 



pn ^(!) pn J 
In r j " 



2ft + A? -A? ' ~" u " n A 4- 



' p« Zi. _ K^>— P n ' 
V/3 ; i 2 " ft 4^ 1 2(l + <) 



_ jr(o)£i 

n ft. 

+ ^n 0) E 



E 



£j Uj y, P? ln ] PI l + d r j- d] 

" - A™ 1 + d? + 4^ 4ft + A™ - A" 1 + d? 



^ 2ft + A 

pn ~(2)pn 



+*n 0) E 



4ft + Af-A™ 1 + d^ 

pn -(1) pn ~(1) pn 



- (2ft + A?-A?)(4ft + A?-A£) i + < 



1 - 



^(Q 2 ) 
« s (p 2 )J 
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(2.34) 



where = The information on the parameters is given in Appendix A. See also 

Ref. [25]. 

Finally, since q^(t) = q^°\t) + <? 7{1) (t) + q l{2) {t), and from (2.5), the moments of the 
flavour-singlet quark, gluon and flavour-nonsinglet quark distributions in the virtual photon 
are given, respectively, by 

qg(n,Q 2 ,P 2 ) = (1, 1) component of the row vector q 7 (t) , (2.35) 
G 7 (n, Q 2 , P 2 ) = (1, 2) component of the row vector q 7 (t) , (2.36) 
qJf S (n,Q 2 ,P 2 ) = (1,3) component of the row vector q 7 (t) . (2.37) 

For n = 2, one of the eigenvalues, A™ =2 , in Eq.(2.30) vanishes and we have <i™ =2 = 0. This 
is due to the fact that the corresponding operator is the hadronic energy-momentum tensor 
and is, therefore, conserved with a null anomalous dimension [5] . The second moments of the 
singlet quark and gluon distributions, qg(n = 2,Q 2 ,P 2 ) and G 1 (n = 2,Q 2 ,P 2 ), have terms 
which are proportional to -^=? and thus diverge. However, we see from (2.33) and (2.34) 



?2 ' 



which vanishes. In 



that these terms are multiplied by a factor 1 — (a s (Q 2 ) / a s ( 
the end, the second moments qg(n = 2, Q 2 , P 2 ) and G 7 (n = 2, Q 2 , P 2 ) remain finite [14]. 



III. RENORMALISATION SCHEME DEPENDENCE 

The structure functions of the photon (nucleon) are expressed as convolutions of coeffi- 
cient functions and parton distributions of the target photon (nucleon). But it is well known 
that these coefficient functions and parton distributions are by themselves renormalisation- 
scheme dependent. There are two kinds of renormalisation-scheme dependence: (i) One is 
the dependence on the renormalisation-prescription (RP) chosen in defining the QCD cou- 
pling constant a s . (ii) The other is the so-called factorisation-scheme (FS) dependence. The 
coefficient functions and parton distributions (equivalently, the anomalous dimensions and 
photon matrix elements of operators) are dependent on the FS adopted for defining these 
quantities. Of course, the physically measurable quantities such as structure functions are 
independent of the choice of the RP for a s and also of the choice of the FS. 

In this section we will show first that the parton distributions of the virtual photon up 
to the NNLO, which were obtained in Eqs.(2.32)-(2.34), are independent of the RP chosen 
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to define a s . Then we consider the parton distributions in the virtual photon target in the 
two factorisation schemes, MS and D1S 7 schemes. 



A. Independence of the renormalisation prescription for a s 

We show that the parton distributions of the virtual photon up to the NNLO are in- 
dependent of the choice of the RS in defining the QCD coupling constant a s . The beta 
function /3(g) is expanded in powers of g 2 up to the three-loop level as 



;Po 



701 



9' 



16tt 2 ^ u (16vr 2 ) 2 11 (16tt 2 ) 3 ' 
Then the QCD running coupling constant a s (Q 2 ) is expressed as [34], 



(3.1) 



P 



a s (Q 2 ) _ J 

4tt ~ {3 L (P L) 2 Po 



-lnL + 



Pi 

(PoL) 3 {(3 



InL- - 



+ 



P0P2 
PI 



o 4r 



i 

L 4 



(3.2) 

where L = ln(Q 2 /A 2 ). It is known that the first two coefficients (3 and Pi in (3.1) are 
renormalisation prescription independent but the coefficient p 2 is not [35]. 

Suppose that p 2 is obtained in one scheme, for example, in the momentum-space sub- 
traction (MOM) scheme [36], and let 5p 2 be a difference between p 2 \one scheme and the one 
calculated in the modified minimal subtraction (MS) scheme [37], 



5P 2 = P. 



2 1 one scheme 



-P. 



2lMS 



(3.3) 



Then we find from (3.2) that the change in the renormalisation prescription for a s (in other 
words, the change Sp 2 ) has an effect on the running coupling constant a s (Q 2 ) as 

'a s (Q 2 ) 



4tt 



(P 



)= 


'« S (Q 2 )1 




47T 



-SP 2 + 0(a 4 s ) , (3.4) 
Po 



so that the effect on the order appears in the higher orders of +2 . 

The change in the renormalisation prescription for a s leads to the change of the coefficient 
functions, the anomalous dimensions and photon matrix elements of hadronic operators. 
We see from (2.32)-(2.34) that the parton distributions in the photon up to the NNLO 
are expressed in terms of the anomalous dimensions of the hadronic operators calculated 
up to the order a 3 (the three-loop level), the mixing anomalous dimensions between the 
photon and hadronic operators up to the order aa 2 (the three-loop level), and photon 
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matrix elements up to the order aa s (the two-loop level). Therefore, only relevant is the 
three-loop anomalous dimension matrix in the hadronic sector. To see this, we use (3.4) and 
we obtain, 

>./3 



(a/) = ^v'W + 



^L^(l) + ^(2) + 



Utt Utt/ A, / ln 



+ 



a 



(An) 



7„ w + 



(4tt) ; 



«_tf(2) 



7„ w + 



;(°) 



4vr /n (4tt) 



« 2 ^/n :i 



~(1) + _Jf-~(2) + . . . 

2 "' (47r) 3 



Thus we get 



In 



^(0) -'(1) = -(1) 
/n i In In i 



^(2)=9'(2)_d2) = ^(0)^2 
" /n In in /n zj 

Po 



(3.5) 



(3.6) 



Let us write the difference of the parton distributions in the photon between one scheme 



chosen to define a s and MS scheme as 



i = 0,l,2, 



(3.7) 



where i = 0,1,2 denote the LO, NLO and NNLO expressions, respectively. Now we find the 
difference in the LO expression q^°\t) given in (2.32) is 



a 



8np l 



An 



d': 



l - 



a s (Q 2 
a s (P 2 



1+dV 



K (0) y- pn _ 



1 



a s {Q 2 ) 5fc 
4tt A) 



q 8 (Q 2 ) ' 
K(^ 2 ). 

« s (p 2 ; 
mq 2 ) 

a s (P 2 ) 



-u- 



where we have used the following formulae: 

An x a s (Q 2 ) 1 



a s (Q 2 
a s (Q 2 ) 



An A) 
a s (Q 2 ) 



a s (<5 



2 U2 



, (3.8) 



(3.9) 
(3.10) 



An J \ An J ) j3 

On the other hand, we see that (3 2 and ^n^ do not appear in the NLO expression q 7 ^(i) 
given in (2.33). Also we already know the effect of the change in (3 2 on the order appears 
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in the higher orders of af +2 . Thus, as far as the analysis for the parton distributions up to 
the NNLO is concerned, we conclude 

6qiW(t)=0. (3.11) 

Finally, we notice that ft and 7^ appear in the NNLO expression of q^ 2 \t) given in (2.34). 
Therefore, we get 



a 



8vrft 



47T 



I ft ^ i-d?V 2) 



(0) ^2 \^ pn a i 

" ft y j 2(1 + 



Pf ffffliy 1 
4ft + A™ - A™ 1 -d n % 



1 - 



1 



^ 4ft + A™ -A? l + < 



«,(g 2 
« S (P 2 ) 



a s {Q 2 ) 
a s (P 2 ) 



d?+l 



Now Eq.(3.6) and the properties of the projection operators in (2.31) give 



P™ SffiP? 



E 



pn g-Wpr. 



4ft + A" - A™ Y 4/90 + X *~ A i 



-p; 



and we obtain 



a 



87Tft 



« s (Q 2 )i5ft 



47T 



4ft ft 

^(Q 2 y 

a s (P2)_ 

« s (g 2 ) 



« s (p 2 ) 



which exactly cancels the change 5qr 7 ^(t) given in (3.8). Thus we find 

5q' (0 \t) + c^ 1 ^) + 5q^ 2 \t) = . 



cG l +l 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



The parton distributions in the photon up to the NNLO are, indeed, independent of the 
renormalisation-prescription adopted to define the QCD coupling constant a s . 

B. Factorisation schemes 



The structure functions are expressed as convolutions of parton distributions and coeffi- 
cient functions. The Mellin moments of the virtual photon structure function ^F^x, Q 2 , P 2 
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is expressed as 

P dxx n - x - P 2 7 (x,g 2 ,P 2 ) = P 2 7 (n,g 2 ,P 2 ) 
Jo x 

= g>,g 2 ,P 2 )-C 2 (n,g 2 ) , (3.16) 

where q y (n, Q 2 , P 2 ) is a four-component row vector 

g>,Q 2 ,P 2 ) = (g 7 (n,g 2 ,P 2 ),^(n,g 2 ,P 2 ),g^(n,g 2 ,P 2 ),r 7 ) , (3.17) 

with = 1, the moment of the photon distribution function (see Eq.(2.2)), being added to 
the row vector q 1 of (2.5), and C 2 (n, Q 2 ) is a four component column vector 

C 2 (n,Q 2 ) = (C|(n,g 2 ), C 2 G (n,g 2 ), Cf 5 (n,g 2 ), C 2 >,g 2 )) T (3.18) 
= (C 2 (n,g 2 ), C 7 (n,g 2 )) T , 

where the hadronic coefficient functions C 2 = (Cf , C^, C-f 5 ) are made up of the flavour- 
singlet quark, gluon and nonsinglet quark, and C 2 is the photonic coefficient functions. 

Since P 2 7 is a physical quantity, its Mellin moments F%(n, Q 2 , P 2 ) in (3.16) is unique and 
FS-independent. But there remains a freedom in the factorisation of P 2 7 into q 1 and C 2 . 
Given the formula (3.16), we can always redefine q 1 and C 2 as follows: 

q\n,Q 2 ,P 2 ) -> q>(n,g 2 ,P 2 )| a = ^(n,g 2 ,P 2 ) Z a (n,g 2 ) , (3.19) 
C 2 (n,g 2 ) - CafaQ^la^Z-^Q 2 ) C 2 (n,g 2 ) , (3.20) 

where q(n,Q 2 , P 2 )\ a and C 2 (n, g 2 )| a correspond to the quantities in a new factorisation 
scheme-a. It is noted that coefficient functions and anomalous dimensions are closely con- 
nected under factorisation. In the following we will study the parton distributions in the 
virtual photon target up to the NNLO in two factorisation schemes, namely, MS and DIS 7 
schemes. 

1. The MS scheme 

This is the only scheme in which the relevant quantities, namely, the j3 function pa- 
rameters up to three-loop level [32], anomalous dimensions up to three- loop level [22-24] 
and photon matrix elements up to two-loop level [38], were actually calculated. We insert 
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them into the formulae given by Eqs.(2.32)-(2.34) and obtain the moments of the parton 
distributions predicted by MS scheme. 

There is another way, a simpler one indeed, to obtain <? 7 (n, Q 2 , P 2 )|ms u P to the NNLO, 
once we know the expression for the moment sum rule of F^(x, Q 2 , P 2 ) obtained up to the 
NNLO and all the quantities in the expression are the ones calculated in MS scheme. We 
will show it in Appendix B. 

2. The DIS 7 scheme 

An interesting factorisation scheme, which is called DIS 7 , was introduced some time ago 
into the NLO analysis of the unpolarised real photon structure function F^x, Q 2 ). Gliick, 
Reya and Vogt [39] observed that, in MS scheme, the ln(l — x) term in the one-loop photonic 
coefficient function Cl{x) for K7, which becomes negative and divergent for x — > 1, drives the 
'pointlike' part of P 7 to large negative values as x — > 1, leading to a strong difference between 
the LO and the NLO results for P 2 7 p 0intlike in the large- a; region. They introduced the DIS 7 
scheme in which the photonic coefficient function C 7 , i.e., the direct-photon contribution 
to is absorbed into the photonic quark distributions. A similar situation occurs in the 
polarised case, and the DIS 7 scheme was applied to the NLO analysis for the spin-dependent 
structure function gJ(x,Q 2 ) of the real photon target [10]. 

The transformation rule from MS scheme to DIS 7 scheme was derived up to the NNLO 
in Ref . [8] . The moments of the parton distributions in DIS 7 scheme are obtained as follows. 
In this scheme, the hadronic coefficient functions are the same as their counterparts in MS 
scheme, but the photonic coefficient function is absorbed into the quark distributions and 
thus set to zero, 

C 2 (n,Q 2 )| DIS7 = C 2 (n,Q 2 )\ m , C 2 >, Q 2 ) | DI s 7 = 0. (3.21) 

Then Eq.(3.16) gives 

F2(n,Q 2 ,P 2 ) = ^(n,Q 2 ,P 2 )| DI s 7 -C 2 (n,g 2 )| DIS7 

= q\n,Q 2 ,P 2 )\ Dl s 7 -C 2 (n,Q 2 )\ m • (3-22) 

On the other hand, F^n, Q 2 , P 2 ) is expressed in MS scheme as 

P 2 >, Q 2 , P 2 ) = qr>, Q 2 , P 2 ) |ms • C 2 (n, Q 2 ) \ m + C^n, Q 2 ) \ m . (3.23) 
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The expansion is made for q J (n, Q 2 , P 2 )|dis 7 in terms of the LO, NLO and NNLO distri- 
butions as 

qr(n, Q 2 , P 2 )|ms 7 = ^ (0) + ^ (1) |dis 7 + ^ (2) |dis 7 + • • • , (3.24) 

where the LO qZ^ is FS-independent. Denoting the difference of <jfn^|Dis 7 (I = 1)2) from 
MS scheme predictions as 5qZ®\ms y i we write 

<£ (0 |dis 7 = q^lm + ^ (0 lms 7 , 1 = 1,2. (3.25) 

Also the hadronic and photonic coefficient functions, C 2 (n,Q 2 )\^ and (n, Q 2 )\ms, are 
expanded in powers of a s (Q 2 ) up to the NNLO as 

r (n o 2 ) I— - r {0) 4. °^l r {1) I— + a ^^ 2 ) r (2) I— + a 9G\ 

*~<2{n,L t > )\ MS - 2n H — — ^ 2 ,nlMS+ (4.7r) 2 2 > ra ' MS ' (o.^o) 



Cl(n,Q 2 )\m = £3n / (e 4 ){cg ) |MS + ^? c 2 ) ^ + ---} ' ^ 

where (e 4 ) = X^i e ^/ n / anc ^ i s FS-independent. Now putting (3.24)-(3.26) into the 
r.h.s. of (3.22) and comparing the result with (3.23) and (3.27), the following relations are 
obtained, 

^in,(e')^\ m = SqZ<% ISl ■ C$ , (3.28) 

SffitW&m = «4>W ■ Cg! + -^"U, ■ ^GgU . (3.29) 

The LO and the NLO C^l^g are written as 

Cg = ((e 2 ),0,l) T , (3.30) 
- (VV 9(1) I— fc 2 V G(1) l— r q{1) \— V (331) 

°2,nlMS — I \ e / C 2,n IMS' \ e / C 2,n IMS' C 2,n I MS J ' ^O.Oij 

where (e 2 ) = ^ f e 2 /n/. Now dividing the quark-charge factor (e 4 ) into two parts, the flavour 
singlet and nonsinglet parts, as 

(e 4 ) = (e 2 )( e 2 ) + ((e 4 )-(e 2 ) 2 ) , (3.32) 

the quark sectors of the difference 5qZ^ |dis 7 at the NLO are given by 

^%s 7 = ^n f (e 2 )c^\ m , (3.33) 
>is 7 = ^3n f ((e 4 ) - (e 2 ) 2 ) c^\ m , (3.34) 
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while we cannot tell anything about 5G^^|dis 7 , since — • 

Actually the DIS 7 scheme was introduced from the very first so that in this scheme the 
photonic coefficient function , i.e., the direct-photon contribution to may be absorbed 
into the quark distributions but not into gluon distribution. Thus we set 5G^|dis 7 = in all 
orders. In other words, we have 

G>, Q 2 ,P 2 ) | DIS7 = G^(n,Q 2 ,P 2 )\ m . (3.35) 

At the NNLO, one obtains from Eqs.(3.29)-(3.34) 

(fo 7(2) U - aa ^ Q2 \ r,{r 2 )(r^ 2) \ r 7(1) I— r q{1) I— f3 36^ 

°QS, nlDIS 7 — (4ti-)2 6H f\ e /^ C 2,n I MS C 2,n I MS C 2,n \MSj > [6.60) 

X„7(2) | _ aa a (Q 2 ) o 4X _ / 2\2\ (7(2). r f(i)i_ |_ ^ <^ 

^ATS, nlDIS 7 — 6 f \\ e I \ e I > ^ C 2,n I MS C 2,n I MS C 2,n ImsJ ■ 1,0.0^ 

These expressions for <5<?s^|dis 7 and „|dis 7 were first derived in Ref. [8]. (See (4.19) 

of Ref. [8]). The parameters c^l^g and \ms are gi ven i n Ref. [5] and c^Ims m Ref. [8]. 
They are also enumerated in Sec. Ill of Ref. [25]. With the knowledge of <5(?s^|dis 7 and 
5q]fg „|dis 7 (I = 1, 2), the parton distributions in DIS 7 scheme <? 7 (n, Q 2 , P 2 )|dis 7 is obtained 
from Eq.(3.25) up to the NNLO. 

Again there is another way to get 9 7 (n, Q 2 , P 2 )|dis 7 up to the NNLO, once we know the 
expression for the moment sum rule of F^ix, Q 2 , P 2 ) up to the NNLO and all the quantities 
in the expression are the ones calculated in MS scheme. It will be shown in Appendix B. 



IV. BEHAVIOURS OF PARTON DISTRIBUTIONS NEAR x = 1 



The behaviours of parton distributions near x = 1 are governed by the large-n limit of 
those moments. In the leading order, parton distributions are factorisation-scheme indepen- 
dent. For large n, the moments of the flavour singlet and nonsinglet LO quark distributions, 
qp°l and q]^ n , both behave as l/(n In n), while the LO gluon distribution behaves 
as l/(n In n) 2 . Thus, in x space, the LO parton distributions in the virtual photon vanish 
for x — > 1 as 
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The behaviours of qJj S (x, Q 2 , P 2 ) for x — ► 1, both in LO, NLO and NNLO are found to 
be always given by the corresponding expressions for qg(x, Q 2 , P 2 ) with replacement of the 
charge factor (e 2 ) with ((e 4 ) — (e 2 ) 2 ). 

In MS scheme, the moments of the NLO parton distributions are written in large-n limit 

as 



^Ims — ^3n /( e 2 )2^, (4.3) 

^ (1) Ims — £3n /( e 2 )l. (4.4) 

The leading contribution to ?s^|ms f° r l &r g e n comes from the (1,1) component of the 
term (Kn^P™ -^-) of Eq.(2.33), which reduces to K^' n -^-. For large n, K^' n behaves as 
(3n/(e 2 ))y^ !i while d n _ as yln n. Then, in x space, we have near x — 1 

qf\x,Q 2 ,P 2 )\ m » ^n f (e 2 )2 [-ln(l - x)] , (4.5) 
G^(x,Q 2 ,P 2 )\ m » ^Sn^e 2 ) [-lux] . (4.6) 

The NLO quark distributions, both q^ 1 \x, Q 2 , -P 2 )|ms an d Qn's ( x > Q 2 > -P 2 )Ims> positively 
diverge as [— ln(l — x)] for x — > 1, while the NLO gluon distribution, GT^^x, Q 2 , -P 2 )Ims> 
vanishes as [—In x]. 

On the other hand, in DIS 7 scheme, the moments of the NLO flavour- singlet quark 
distribution (?5^|dis 7 is expressed in large-n limit as (see Eq.(3.25)), 

„7(1)| _ 7(1) I 1 X«7(l) I 

nlDIS 7 - fe,n I MS + °<ls, nlDIS 7 



47r" 



j 7 
In n 



n J 



Sn^e 2 ) -2 — - , (4.7) 



since c^Ims i n Eq.(3.33) behaves as (—4 In n)/n for large n. Thus we have for large x, 

qf\x,Q 2 ,P 2 )\ D1Si » ^3n / (e 2 )2 ln(l - x) . (4.8) 

The NLO distribution in DIS 7 scheme, q^\x, Q 2 , -P 2 )|dis 7? negatively diverges as x — > 1. 
This is due to the fact that the NLO photonic coefficient function (the inverse Mellin 

transform of c^), which in MS becomes negative and divergent for x — > 1, is absorbed into 
the quark distributions in DIS 7 scheme [39]. 

The moments of the NNLO parton distributions in MS scheme behave for large n as 

q ^\ m _^ ^^S!l3n f (e 2 ) 8 -^ , (4.9) 
ys,n ims 4vr 4tt /W 9 n ' K J 
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The leading contribution to q^nlus f° r l &r g e n comes from the (1,1) component of the term 



(-K {2) P™ j^r) in Eq.(2.34), which reduces to K ( s f hn ±r. The large-n behaviour of Kg 



(2),n 1 



-(2),n 



is (3n f (e 2 )) 



2\\ 256 (Inn) 4 



27 n 



. On the other hand, the large-n behaviours of 

80 (In n) 



.7(2) 



,7(1) I 



,n IMS' C 2,n' I MS aI1< ^ 



c 2,n I Ms are gi ven by — y ra , — and |(ln n) 2 , respectively, and thus we see from 

Eq.(3.36), 

(4.11) 



7(2) | 
5, nlDIS 7 



aa s , , 2 \A 16 ( ln 



for large n. So we find in DIS 7 scheme, 



^' (2) |dis 7 



9s,n I MS + 6 IS, nlDIS 7 

« ^(^ 2 ) 3n/ ( e 2\ 8( ln «) 3 



7(2), 



An An 1 v ' 3 n 
In x space, therefore, the NNLO parton distributions near x = 1 are 

,7(2), 



qf\x,Q 2 ,P 2 )\ mSj 



G« 2 \x,Q 2 ,P 2 )\ m 



a 


a s (Q 2 ) 


Air 


An 


a 


a s (Q 2 ) 


An 


An 


a 


a s (Q 2 ) 


An 


An 



3n f (t 



9 

8 

3 
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3n f {e 2 )- -\n-\l-x) 



-ln 3 (l - x) 



9 



3nf{e 2 ) — — lnx ln 2 (l — x) 



(4.12) 

(4.13) 
(4.14) 

(4.15) 



It is noted that NNLO quark distributions in both MS and DIS 7 schemes diverge at x = 1 as 



[— ln 3 (l— x)] and, furthermore, that q^ 2 \x, Q 2 , -P 2 )|dis 7 rises sharper than q^^x, Q 2 , P 



7(2), 



IMS 



as x — > 1. This sharp rise of q 1 g 2 \x, Q 2 , P 2 )|dis 7 as [— ln 3 (l — x)] near x = 1 is an unexpected 



result. Indeed, both the NLO and NNLO photonic coefficient functions c^ l \x) and c], [Z> (x) 
in MS scheme negatively diverge as x — > 1 (see above for the large-n behaviours of c^\mE 
and c^Ims). Then, the first thought is that absorbing these photonic coefficient functions 
into the quark distributions would make the DIS 7 quark distributions behave milder than 
those in MS scheme or negatively diverge as x — > I. The NLO quark distributions in DIS 7 
scheme work fine but not the NNLO quark distributions. This is due to the contribution of 
the second term (-c^Ims 4$\m) to <^%s 7 in Eq.(3.36). 



,7(2), 



V. NUMERICAL ANALYSIS 



The parton distributions in the virtual photon are recovered from their moments by the 
inverse Mellin transformation. In Fig. 2 we plot the parton distributions in MS scheme in 
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3GeV A=200MeV (MS) 




n,=4 Q=100GeV z P^Ge^ A=200MeV (MS) 




(a) (b) 

FIG. 2: Parton distributions in the photon in MS scheme for rif = 4, Q 2 = lOOGeV 2 , P 2 = 3GeV 2 
with A = 0.2GeV: (a) xq],(x,Q 2 , P 2 )\ m and (b) xG"< {x , Q 2 , P 2 )\ m . 



units of (3nf(e 2 )a/n)\n(Q 2 /P 2 ): (a) the singlet quark distribution xqg(x,Q 2 , P 
(b) the gluon distribution xG^(x,Q 2 ,P 2 )\ m . We have taken n f = 4, Q 2 = lOOGeV 2 , P 2 = 
3GeV 2 , and the QCD scale parameter A = 0.2GeV. We see that both the (LO+NLO) and 
(LO+NLO+NNLO) curves show the similar behaviours in almost the whole x region, which 
means that the NNLO contribution is small. The behaviours of these two curves, however, 
are quite different from the LO curve. They lie below the LO curve for 0.2 < x < 0.8, but 
diverge as x — ► 1. Compared with the quark distribution, the gluon distribution xG 1 \y^ is 
very small in absolute value except in the small- x region. Concerning the nonsinglet quark 
distribution xq^six-, Q 2 -, P 2 )\uEi we nnc ^ that wnen we take into account the charge factors, 
such as (e 2 ) = Ylii e V n f an d (e 4 ) — ^2i e f/ n f, it falls on the singlet quark distribution in 
almost the whole x region; namely the two "normalised" distributions xq^ = xq^/ie 2 ) and 
xq^ s = xqNs/(( eA ) ~ ( e2 ) 2 ) mostly overlap except at the very small x region. This situation 
is the same in both MS and DIS 7 schemes. The rise of the singlet quark distribution near 
x = is related to the gluon distribution which grows rapidly as x — ► 0. 

The parton distributions in DIS 7 scheme were analysed up to the NNLO. In Fig. 3 (a) we 
plot the singlet quark distribution xq$(x, Q 2 , P 2 )|dis 7 m units of (3n/(e 2 )a/7r)ln((5 2 /P 2 ). 
Again we chose n f = 4, Q 2 = lOOGeV 2 , P 2 = 3GeV 2 , and A = 0.2GeV. The three curves 
(LO, LO+NLO, LO+NLO+NNLO) rather overlap below x = 0.6. Absorbing the photonic 
coefficient function Cl into the quark distributions in DIS 7 scheme has an effect on their 
large-x behaviours: Unlike the MS scheme, the (LO+NLO) curve goes under the LO curve at 
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n,=4 Q 2 =100GeV 2 P 2 =3GeV 2 A=200MeV (DIS V ) n,=4 Q 2 =100GeV 2 P 2 =3GeV 2 A=200MeV (DIS T ) 




(a) (b) 

FIG. 3: (a) Singlet quark distribution xq^(x, Q 2 , P 2 )|dis 7 i n DIS 7 scheme for nj = 4, Q 2 = 
lOOGeV 2 , P 2 = 3GeV 2 with A = 0.2GeV. (b) LO, NLO, and NNLO contributions. 

x ~ 0.6 and the difference between the two grows as x — > 1. Adding the NNLO contribution 
makes the difference bigger at large x except near x — 1. At very close to x — 1 the 
(LO+NLO+NNLO) curve shows a sudden surge. In order to see the details, we plot in 
Fig.3(b), the NLO and NNLO contributions in DIS 7 scheme. We observe that the NLO 
contribution is large and negative for x > 0.8, while the NNLO contribution remains to 
be very small until very close to x = 1 and then blows up. The behaviours near x = 1 of 
g^ (1) (x,(5 2 ,P 2 )| D is 7 and q^ix, Q 2 , -P 2 )|dis 7 were discussed in Sec. IV. 

Finally, the gluon distribution xG^(x, Q 2 , -P 2 )|dis 7 is the same as xG 7 (x, Q 2 , -P 2 )|ms- 

VI. CONCLUSIONS 

We have analysed the parton distributions in the virtual photon target which are predicted 
entirely up to the NNLO in perturbative QCD. Parton distributions are dependent on the 
scheme which is employed to factorise structure functions into parton distributions and 
coefficient functions. The virtual photon target serves as a good testing ground for examining 
the behaviours of the parton distributions and their factorisation-scheme dependences. We 
have studied the quark and gluon distributions in two factorisation schemes, namely, MS 
and DIS 7 schemes. 

We see from Figs.2(a) and 3(a) that (LO+NLO) and (LO+NLO+NNLO) curves for the 
quark distribution xq^(x, Q 2 , P 2 ) show quite different behaviours in two schemes, especially 
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in the large-x region. From the viewpoint of "perturbative stability", the DIS 7 scheme 
gives a more appropriate behaviour for xq1(x,Q 2 , P 2 ) than MS. The gluon distribution 
xG"'(x, Q 2 , P 2 ) is the same in both schemes and is predicted to be very small in absolute 
value except in the small-x region. Finally, we observe that the (LO+NLO+NNLO) curve 
for xq$(x, Q 2 , P 2 )|dis 7 shows a sudden surge at very close to x — 1. Although the NLO 
contribution, |dis 7 , negatively diverges as ln(l — x) for x — > 1, the NNLO (?s^|dis 7 
diverges positively as [— ln 3 (l— x)\. This may hint a necessity of considering the resummation 
for parton distributions and also for the photon structure function F^(x : Q 2 , P 2 ) itself near 
x = 1 [40]. 
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APPENDIX A: PARAMETERS IN q^ 0) (t), g 7(1) (i) AND q^ 2) (t) 



We give here some important information on the parameters which appear in q^°\t), 
<? 7(1) (t) and <? 7(2) (t) in Eqs.(2.32)-(2.34). They are all calculated in MS scheme. We introduce 
the following quark charge factors: 

"V n n f 



(e 2 



fs = 3n f (e 2 } , f NS = 3n / ((e 4 ) - (e 2 ) 2 ) . (Al) 

1. Initial conditions for parton distributions: A n ^ and 

The elements of the row vector A$ are expressed as = (^An^\ An^\ An S ^ j , with 
1 — 1,2. Then we have at one-loop level, 

%V = f s HV(n), i G(1) = 0, A N n ^ = f NS H^\n), (A2) 

where the explicit expression of H^\n) is given in Eq.(3.30) of Ref.[25]. The elements at 
two-loop level are 

%W = f 3 HV(n), A°V = f 3 HS\n), A^ = f NS Hf\n) , (A3) 
where Hf\n) and H^\n) are given in Eqs.(3.39) and (3.42) of Ref.[25], respectively. 

2. Anomalous dimensions 

The one-, two- and three-loop anomalous dimensions for the hadronic sector, 7™ \ In ^ 
and 7r t 2 ' ) [22, 23], respectively, are already in literature. The eigenvalues A™ of the one- 
loop anomalous dimension matrix 7^ and the corresponding projection operators P™ which 
appeared in Eqs. (2.30)-(2.31) are given in Ref.[5]. Also the one- and two-loop photonic 
anomalous dimensions, iff and K n l \ are already known [5, 33, 39]. But concerning the 
three-loop anomalous dimensions K^' n , K^' n and K^ n , the exact expressions have not 
been in literature yet, but approximate ones were given [24]. It is remarked there that the 
precision is within about 0.1% or less. The approximate expressions for Kf ' n , K^ n and 
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K®' n are 

4 2) '" « 4 2) ap pr ox = ~fs 2{i?— (n) + £ P s»} , 

4 2) ' n « 4 2) r ppr ox = -fs (n) , (A4) 

Kg* « ^' n approx = -/at5 2i?— (n) , 

where the explicit expressions of £'^ > P rox (n), E ps ^(n), and i?^ prox (n) are given, respectively, 
in Eqs.(B2), (B4) and (B3) of Ref.[25]. 

APPENDIX B: ANOTHER WAY TO FIND PARTON DISTRIBUTIONS IN MS 
AND DIS 7 SCHEME 

Once we know the expression for the moment sum rule of F^x, Q 2 , P 2 ) obtained up to 
the NNLO corrections [25] and all the quantities in the expression are the ones calculated in 
MS scheme, then there is an easy way to find the parton distributions in the virtual photon 
in both MS and DIS 7 schemes up to the NNLO. In the following the equation numbers 
correspond to those in Ref. 



(i) Parton distributions up to the NNLO in MS scheme 

In the moment sum rule of F 2 7 (x, Q 2 , P 2 ) given by Eqs.(2.29)-(2.37), we set C^l = 
Cfl = and C^n — CjJn^ — 0- In addition, we put C^l = (1, 0, 0) T , then we obtain 
ql{n,Q 2 ,P 2 )\ m . Similarly, when we put eg = (0,1, 0) T and eg = (0,0, 1) T , we 
obtain G y (n, Q 2 , -P 2 )|ms and Q.Ns( n i Q 2 ^ -P 2 )Ims > respectively. 

(ii) Parton distributions up to the NNLO in DIS 7 scheme 

For the gluon distribution, we have G 7 (n, Q 2 , P 2 )|dis 7 = G 7 (n, Q 2 , -P 2 )|ms • 

• qJj S (n, Q 2 , -P 2 )|dis 7 is obtained as follows: 

Winner I „ „ _7 I _„T,(0)| , 7,(1) I , _7,(2) I 

Writing q NS \Dis^ as q N s\Dis-, - q NS |dis 7 + Qns Idis 7 + Qns Idis 7 , 



,7.(0) I 



DIS 7 / 


a 




4-7T 






[*s(Q 2 )\ 



is obtained from the term C% s in Eq.(2.30), which is 



the same with ?jvs^Ims- 



~~ QwP |dis 7 / * s obtained from the sum of the terms ^4^ 5 , £>j^ 5 and C™ in 

Eqs.(2.31)-(2.33), where we put, Cfl = (0,0, 1) T , c£J = and is replaced 
withC^x^. 
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~ Qns Idis 7 / 



MQ 2 ) 

4-7T 



is obtained from the sum of the terms T>^ s , £^ S: 



and g n in Eqs.(2.34)-(2.37), where we put, Cf ] n = (0,0, 1) T , C 
and <) is replaced with (< 2) - < } <f«) x 



2,n 



n r (2) 







<?s(n, Q 2 , -P 2 )|dis 7 is obtained as follows: 
Writing g^| D is 7 as g^Dis, 



DIS-, 



7,(1) 

Qs 



DIS-, 



DIS-, 



<& m < 



DIS 7 / 


a 




4-7T 


_87rA,_ 




.MQ 2 ). 



which is the same as q$'^ 
where we put, C^l = (1, 0, 0) T 



is obtained from the terms (c%+C n _)j /(e 2 ) in Eq.(2.30), 
\ms- Or we can obtain from the sum of the terms J2i 



DIS. 



V 



is obtained from the sum of the terms J2i A™, J2i an d ^™ i n 



Eqs.(2.31)-(2.33), where we put, Cf] = (1,0, 0) T , c£J = and C$ is replaced 



with .C™ x gl. 



7,(2) 
55 



'2,n 
DIS 7 / 



8tt/3 



MQ 2 ) 

4-7T 



is obtained from the sum of the terms V*.Z>J\ 
Ei^T and G n in Eqs.(2.34)-(2.37), where we put, C$ = (1,0, 0) T , C$ = 0, 
= and q ( „ 2) is replaced with (c#> - C™C™V>) x Note 



NS(1) 
2,n 
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